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We study the interface dynamics of a discrete model previously shown [A. Sanchez, M. J. Bernal, 
and J. M. Riveiro, Phys. Rev. E 50, R2427 (1994)] to quantitatively describe electrochemical deposi- 
tion experiments. The model allows for a finite density of biased random walkers which irreversibly 
stick onto a substrate. There is no surface diffusion. Extensive numerical simulations indicate that 
the interface dynamics is unstable at early times, but asymptotically displays the scaling of the 
Kardar-Parisi-Zhang universality class. During the time interval in which the surface is unstable, 
its power spectrum is anomalous; hence the behaviors at length scales smaller than or comparable 
with the system size are described by different roughness exponents. These results are expected to 
apply to a wide range of electrochemical deposition experiments. 



The field of non equilibrium dynamics of rough sur- 
faces and interfaces has undergone an enormous activity 
during the last decade. Considerable experimental and 
theoretical efforts have been devoted to understanding 
the common general features of the growth of surfaces in 
various seemingly unrelated phenomena ■ By exploit- 
ing the existence of universality and scale invariance akin 
to those present in equilibrium critical phenomena, most 
of the activity has been aimed to organizing non equilib- 
rium surface growth processes into universality classes. 
In this context, surfaces are described by the global width 
(or roughness) W(L, t), the rms fluctuations of the height 
variable h(x,t) giving the location of the surface at posi- 
tion x over a substrate of lateral size L at time i, around 
its mean value hr,{t) = x ^2x h(x, t): 

W 2 {L,t) = ±{Y^{h{x,t)~h L {t)f), (1) 

X 

where angular brackets stand for noise average. 

In many cases it is observed that the width initially 
grows as W(L,t) ~ t 13 for t <C L z , saturating at a size 
dependent value W(L) ~ L a for t 3> L z . The rough- 
ness exponent a, the dynamic exponent z and their ratio 
(3 = a/z are usually taken to identify the universality 
class the growth process considered belongs to. Under 
the assumption of scale invariance (self- affinity), analo- 
gous behavior is expected for the rms fluctuations of the 
height within a box of lateral size I <C L, or local width 

w 2 (l,t) = ^(£(Kx,t)-h l (t)f), (2) 

X 

namely 

»(M)~(j„ if t>>P . (3) 



The scaling behavior (|^) holds, e.g., in the Kardar-Parisi- 
Zhang (KPZ) equation ||, which generically describes 
surfaces growing irreversibly in the absence of specific 
conservation laws However, in recent studies ||,||] of 
surface growth in the context of Molecular Beam Epitaxy 
(MBE) and related vapor deposition techniques, it has 
been found that many relevant discrete models and con- 
tinuum equations do not obey the local scaling (||), but 
rather for intermediate times l z <C t <C L z it is found 
that w(l,t) ~ l a >°" if"- ""''/ 2 , with a ioc / a; hence 
the behaviors of the surface at small and large length 
scales differ. This phenomenon has been termed anoma- 
lous scaling, and is not yet understood at a fundamental 
level (see HQ for overviews). In principle, anomalous 
scaling was somehow associated with super-roughening 
(i.e., the existence of exponent values a > 1). However, 
recent works have shown that it can appear as an inde- 
pendent phenomenon |Jj (henceforth referred to as in- 
trinsic anomalous scaling) caused by anomalous behav- 
ior of the surface power spectrum |^|;[l0| (to be defined 
below). Up to now, the only cases in which anomalous 
scaling has been identified are related to models with a 
conserved current along the surface, as those for MBE 
growth 1^,0], or with some type of disorder J8|-[l^|. We 
are aware of no previous report on the occurrence of in- 
trinsic anomalous scaling in models with time dependent 
noise and no constraint or conservation law. Moreover, 
an obviously important open issue is to clarify the exper- 
imental conditions under which anomalous scaling can be 
observable. 



In this Rapid Communication we study a surface 
growth model, called multiparticle biased diffusion lim- 
ited aggregation (MBDLA), introduced in |e|] and 
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shown to describe quantitatively electrochemical depo- 
sition (ECD) experiments Oj. In particular, it allows to 
study systems with a finite density of depositing parti- 
cles of different species. The relaxation rules of MBDLA 
do not yield a conserved current. Specifically, consistent 
with a range of experimental applications [13] , there is no 
surface diffusion. As a consequence, we will show below 
that MBDLA asymptotically displays the scaling behav- 
ior of the KPZ universality class. Nevertheless, there 
does exist a time regime in the evolution of the surface 
during which the effective scaling is intrinsically anoma- 
lous and described by non universal exponents. This time 
regime is associated with the onset of an instability due 
to the Laplacian character of the model. This feature is 
analogous to recent reports jjJI that for a class of model 
discretized equations in which there exists a conserved 
current, anomalous scaling with non universal exponents 
is present during a transient regime in which the sur- 
face growth is unstable. Here we find an analogous phe- 
nomenon for a realistic model without conserved currents 
which belongs to the KPZ universality class. Actually, 
our findings are reminiscent of a discrete non conserved 
erosion model fll6|] also displaying an initial instability 
followed by asymptotic KPZ scaling. The crucial differ- 
ence is that the model in |l6) is described by the noisy 
Kuramoto-Sivashinsky equation, asymptotically equiva- 
lent to KPZ in 1+1 dimensions 0, and — as KPZ itself — 
free of anomalous scaling. 

The algorithm that defines MBDLA begins with 
a number of random walkers (cations) randomly dis- 
tributed with concentration c on a two-dimensional 
square lattice with periodic boundary conditions. The 
lower side of the lattice is chosen to be the cathode. 
The initial condition evolves in time as follows: Every 
time step a walker is chosen and moved to one of its four 
neighboring sites with different probabilities: 0.5 to move 
parallel to the cathode (either left or right); 0.25 + p to 
move down towards the aggregate (the cathode, at time 
t = 0) and 0.25 — p to move up, away from it. The pa- 
rameter p is referred to as the bias and, as shown in [fl3| , 
can be quantitatively related to the electric current in 
the physical system; experimentally, it can be controlled 
by changing the value of the applied electric field. After 
a destination site has been selected, the particle moves 
if that node of the lattice is empty (and we select an- 
other particle if it is not). Once the particle has been 
moved, if the new position has any nearest neighbor site 
belonging to the aggregate, the walker's present position 
is added to the aggregate with probability s; otherwise it 
stays there (and is able to move again) with probability 
1 — s. We term s the sticking probability; it is related to 
the chemical activation energy the cation needs to stick 
to the aggregate. 

Of the two parameters of MBDLA, the bias p and the 
sticking probability s, only p is independent for the ex- 
perimental description, s being simply proportional to it 



|13| . This remark notwithstanding, in this work we will 
consider both p and s as independent parameters, for 
our goal is to make a general analysis of the model be- 
havior beyond the experimentally motivated values. As 
seen in what follows, this will help to obtain a coher- 
ent, complete picture of MBDLA properties. Thus, s 
takes values between and 1, and p ranges from to 
0.25. Finally, let us note that MBDLA produces in gen- 
eral a non compact aggregate which can have voids and 
overhangs. Large values of p produce nearly ballistic 
trayectories, while smaller values drive the particles in 
a Brownian motion fashion. In two limits, the behavior 
of MBDLA approaches other relevant models: ballistic 
deposition (for p = 0.25), well known to belong to the 
KPZ universality class and not to display anomalous scal- 
ing [^J, and multiparticle Diffusion Limited Aggregation 
(MDLA) (for p — 0) , the paradigmatic model of unstable 
Laplacian growth jjj . 

To begin the summary of our results, we will describe 
the phenomenology for p = 0.05 and arbitrary s. Phys- 
ically, p — 0.05 represents not too small applied electric 
fields, surface diffusion being negligible, so that MBDLA 
provides a good description of ECD. We will later dis- 
cuss how changes in p modify the description that fol- 
lows. Simulations were performed in a L — 300 by 1500 
lattice, although other values of L were also considered 
(see below). The dynamics of the model is conveniently 
characterized by the global width, Eq. ([!]). 
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FIG. 1. Global width vs time for p = 0.05, s = 1, c = 0.1. 
The dynamics presents an initial transient during which the 
surface is still uncorrelated and /3 ~ ^. Then there is an 
unstable transient within which (5 ~ 0.7. In this region the 
system presents intrinsic anomalous scaling (see text). Be- 
fore saturation, the KPZ nonlinearity governs the dynamics, 
hence (3 ~ 0.33. The present plot is an average over 50 noise 
realizations. All magnitudes are measured in arbitrary units. 

We are interested in the dynamics of the active zone 
|EJ, hence we define the variable h(x,t) to be the height 
of the topmost particle belonging to the cluster for each 
column. As mentioned above, for times before satura- 
tion due to finite system size (t < L z ), the width is ex- 
pected to scale as W(L,t) ~ t^ . In our case, we find 
that, along the time evolution, (3 takes up to three differ- 
ent values for large enough systems, see Fig. 111. At early 
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times, shot noise dominates the growth, thus (3 ~ 0.5 
as in a random deposition process. As time proceeds, 
the Laplacian character of the model shows up, favoring 
the unstable development of pillars, which in turn leave 
deep grooves behind. This instability is reflected in the 
large value of the growth exponent ((3 ~ 0.59 — 0.73 in- 
creasing with sticking probabilities). As local slopes of 
pillars (or grooves) become larger, they seem to trigger 
further nonlinear effects, whose consequence is the sta- 
bilization of the surface dynamics. Thus, for s > 0.5 
we measure in this regime (3 ~ 0.25, a value close to 
that of the Edwards- Wilkinson universality class [ fL8| . 
For smaller values of the sticking probability, s < 0.5, 
we measure (3 ~ 0.33, the growth exponent characteris- 
tic of the KPZ universality class. Actually simulations 
for larger system sizes (L — 512, 1024) yield (3 ~ 0.33 
before saturation for almost all values of s. Hence the 
asymptotic behavior of the model is described by KPZ 
scaling. This also shows the role of s as a noise reduc- 
tion parameter flTi| , large noise reduction (s — > 0) tuning 
the system closer to the asymptotic scaling regime. The 
value of the roughness exponents are ai oc ~ a ~ 0.5 in 
all cases, showing that there is no anomalous scaling in 
the asymptotic state. Measurements of the mean excess 
velocity v(m) of an interface subject to an average tilt 
to, imposed through helical boundary conditions, confirm 
fl9| the relevance of the KPZ nonlinearity in the effec- 
tive description of our system. Once established that the 
long time scaling behavior of the interface is of the KPZ 
type, we focus on the unstable transient regime, char- 
acterized by two main features, namely a large (3 value 
(hence very rapid growth) and intrinsic anomalous scal- 
ing caused by a non standard form of the power spec- 
trum of the surface, S(k,t) — (h(k,t)h(—k,t)}, where 
h{k, t) = L- 1 ' 2 EJH^ *) - M*)] exp(ifcE). 




FIG. 2. Collapsed power spectrum within the unstable 
transient using the same parameters as in Fig. [jj. The in- 
set shows the time evolution of the power spectrum every 10 6 
time units. Thick lines have slopes given by the exponents 
quoted in Table I. All magnitudes are measured in arbitrary 
units. 

In Fig. U we plot the colapse of the power spectrum 
S(k,t) within the unstable transient, where the inset 



shows S(k,i) as function of k for every 10 6 time units. 
S{k, t) displays a behavior consistent with the scaling 
form 

S(k,t) = k-<> 2a+1 h{kt l / z ), (4a) 

where 

■< u ) = \ u ^ if U « L ( 4b ) 

The exponent 9 — a — on oc measures the difference be- 
tween the global roughness exponent a and the local 
roughness exponent ai oc . The exponents are a ~ 2.15, 
ctioc ~ 0.5, (3 ~ 0.73, and z ~ 2.95 (all exponent val- 
ues given are to within an accuracy less than 5%). It has 
been shown in Ref. Jl(| that a structure factor such as (jij) 
implies intrinsic anomalous scaling for the local width. 
The slopes of the collapsed curves yield exponents con- 
sistent with the ones used to achieve the best collapse. 
Collapses of structure factor data obtained for different s 
values yield similar qualitative results, as shown in Fig. |^, 
and the exponent values summarized in Table I. All the 
exponent values (except notably a/ oc , see also p|fl^0| ) 
depend strongly on model parameters, and we consider 
them as non universal effective values. 




FIG. 3. Collapsed power spectra within the unstable tran- 
sient using the exponent values in Table I, with s — 0.75 and 
s = 0.25. Parameters other than s are as in Fig. |l|. Thick 
lines have slopes given by the exponents quoted in Table I. 
All magnitudes are measured in arbitrary units. 

To provide a general description of our model, we have 
studied bias values different from p — 0.05. For example, 
for p — 0.005 and large values of the sticking probabil- 
ity s, we have observed DLA like patterns as expected 
fill . Consequently the unstable regime lasts longer than 
in the cases studied above. The same is observed for 
larger values of s at low p values. Moreover, decreasing 
(increasing) the ion concentration c increases (decreases) 
the duration of the unstable growth, in accordance with 
the fact that as c — > pure DLA is recovered. Hence 
we can tune the two parameters and the concentration 
to reproduce different experimental behaviors, because 
changes in the parameters alter the relative duration of 
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the different scaling regions in the dynamics of the sys- 
tem. However the asymptotic KPZ scaling still holds at 
long enough times. 

In summary, we have shown the existence of intrinsic 
anomalous scaling for a realistic model, MBDLA, with- 
out conserved currents and hence in the KPZ univer- 
sality class. Anomalous scaling appears within a time 
region of variable duration (depending on parameter val- 
ues) in which the surface evolves in an unstable fashion. 
This type of scaling is due to a non-standard form of the 
power spectrum and features non-universal (parameter 
dependent) values of the effective exponents. We inter- 
pret anomalous scaling to be associated with the unstable 
structures developing in the model dynamics. Thus, for 
large length scales we measure a > 1, associated with the 
steep slopes of pillars. On the other hand, when probing 
short length scales we obtain ai oc ~ 0.5 (see Table I), a 
value numerically close both to the KPZ exponent and 
to the exponent resulting from the convolution of the 
Heaviside function with the interface pi]] (through our 
definition of the height variable h(x,t)). At any rate, in- 
trinsic anomalous scaling seems to be strongly associated 
in our case with the nonlocality of the evolution rules of 
the model which produces the unstable structures of the 
interface. 

To date, anomalous scaling has been mostly observed 
in theoretical models, hence the interest of our results: 
as MBDLA is a correct description of ECD [[u]], our 
findings could be observed in that kind of experiments. 
Interestingly, MBDLA displays instabilities similar to 
those reported in various ECD experiments [ pp| , reinforc- 
ing the connection of the model to actual physical sys- 
tems. Moreover, anomalous scaling may explain the wide 
spread in the a values (a ~ 0.5 — 0.8) reported in various 
ECD experiments p0| . As shown in ||, in the presence 
of anomalous scaling, the fact that all length scales sat- 
urate at the same global saturation time t sat ~ L z poses 
additional difficulties to the evaluation of exponents by 
means of the local width. Specifically, if not performed at 
global saturation, the roughness exponent thus obtained 
is an effective value different from both, a and ai oc . Re- 
search on extensions of the model to more experimental 
parameters is under way, the influence of surface diffusion 
being the center of our future work. 
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P 


Q 


Otloc 


z 


1.00 


0.73 


2.15 


0.50 


2.95 


0.75 


0.67 


2.04 


0.49 


3.03 


0.50 


0.67 


1.92 


0.54 


2.87 


0.25 


0.59 


1.73 


0.50 


2.94 



TABLE I. Summary of relevant exponents as measured 
with S(k,t) and w(l,t) within the unstable transient, for 
p = 0.05, s = 1, c= 0.1. 
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